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We present measurements of temperature and magnetic field dependence of the critical current and 
excess current in a carbon nanotube Josephson quantum dot junction. The junction is fabricated in 
a controlled environment which allows for extraction of the full critical current. The measurements 
are performed in the open quantum dot regime, and fitted to theory with good qualitative agreement. 
We also show how to extract level spacing, level broadening, and charging energy of an open quantum 
dot from a bias spectroscopy plot. 
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Nanoscale Josephson quantum dot junctions are in- 
triguing devices showing several interesting physical phe- 
nomena. Supercurrent, Andreev reflections, quasiparti- 
cle transport, and excess current have all been studied in 
junctions where a nanotube or nanowire constitute the 
quantum dot [H [2 O gHH El [7] . Furthermore, the inter- 
play between these Josephson junction related phenom- 
ena and correlations as the Kondo effect [HI [9l [TOl HTJ [T^ 
and the O-tt transition for more weakly coupled junctions 
has been explored [131 IHl US] • 

In this paper, we present experimental results in the 
strongly coupled regime for a Josephson quantum dot 
junction realized in a carbon-nanotube. Inspired by 
Ref.[16 , we utilize a designed external circuit in order to 
control the phase fluctuations which enables us to infer 
the true magnitude of the critical current, /c, from the 
measurable critical current /switching current, /^n, by a 
fitting procedure [13 . Im can significantly differ from Ic 
as demonstrated previously for nanotube-based Joseph- 
son junctions [3l jH |13]. Here we analyze the magnetic 
field dependence and temperature dependence of both 
the critical current and excess current. 

The devices are fabricated on a degenerately doped sili- 
con wafer with a 0.5 /im layer of Si02. Carbon nanotubes 
are grown from islands of catalyst material and contacted 
by small electrodes of superconducting trilayers of 5nm 
Ti, 60nm Al and 5nm Ti. The superconducting elec- 
trodes are kept small to reduce junction capacitance. 
Each superconducting electrode is contacted by two nor- 
mal metal leads to bonding pads which enables four probe 
measurements. The measurements are performed in a 
^He-^He dilution fridge with a base electron temperature 
of 75 mK. Inside the dashed square in Fig.jlja) we show 
a schematic circuit diagram of the on-chip components of 
the full Josephson junction. The fabrication is similar to 
Ref. iMi. The superconductor-nanotube-superconductor 
junction is represented by a Josephson element (cross), a 
junction resistor and a junction capacitor Cj. The 
Josephson element has a current-phase relation, which 
we in the fitting procedure (see below) assume to be 
/((/)) = /c'sin((/)), with (j) being the phase difference be- 
tween the two superconducting and Ic the critical cur- 




FIG. 1: (color online) (a) Inside the dashed square: Schematic 
circuit of on-chip components of the Josephson junction. Out- 
side dashed square: Four probe voltage bias setup for mea- 
suring junction voltage Vj vs. current /. The full Joseph- 
son junction consists of both the superconductor-nanotube- 
superconductor junction, represented by a Josephson element 
Ic in parallel with a junction capacitor Cj and junction resis- 
tor and on-chip resistors R and capacitances C. (b) Bias 
spectroscopy plot of differential conductance versus source- 
drain, and gate voltage, (c) Schematic of a Fabry-Perot dia- 
mond, (d-f) Schematic transport diagrams at zero-bias reso- 
nance (d), positive-bias resonance (e), and negative-bias res- 
onance (f). 



rent. However, the sinusoidal form of this relation is not 
in general true and this may cause some inaccuracy in 
the determination of the critical current. In Ref. [13 we 
show that the difference between the extracted Ic using 
either sin(^) or the correct functional form of the current- 
phase relation is in fact small and moreover largest near 
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the O-TT transition relevant only for closed dots. We there- 
fore expect the simpler relation also to be a reasonable 
approximation in the case of open dots (which allows 
us to use the Ivanchenko-Zil'berman relation in Eq. (|9|). 
At sub-gap bias voltages Rj accounts for current due to 
multiple Andreev reflections and at higher bias voltages 
it accounts for quasi-particle transport. The capacitance 
between the superconducting electrodes {Cj ~ 5fF), and 
between bonding pads (C ~ 1 pF) is estimated as a par- 
allel plate capacitance through the back gate. We have 
fabricated long thin metal leads with a measured resis- 
tance of ~ 1 kl] between the bonding pads and the su- 
perconducting electrodes which, as will be shown later, 
is crucial for increasing the measurable critical current. 

In Fig.jljb) we show a bias spectroscopy plot of differ- 
ential conductance versus source-drain voltage {Vsd)-, and 
gate voltage iVgate)' Regular conductance oscillations in 
both source-drain and gate voltage is seen due to tun- 
ing of successive energy levels in the dot, with a sepa- 
ration (level spacing) A^, on and off resonance. Bias 
spectroscopy plots with the leads in the normal state 
{B = 150 mT) (not shown) show conductances at the 
resonances ranging from 2 to 3.5e^/h. That high con- 
ductance is only allowed when the degeneracy of each 
energy level is four- fold (spin and orbital), and when 
the broadening of each energy level F = + F^^, where 
the Tg/h {Td/h) is the tunnel rate through the source 
(drain) barrier, is larger than the Coulomb repulsion en- 
ergy for adding an electron to the dot (charging energy) 
Uc = e^jC. This regime (A^ > F > /7c) is often called 
the Fabry-Perot regime [3j ilTj , and the dot is termed an 
open quantum dot. We will now analyze the bias spec- 
trum and extract energy parameters, tunnel couplings, 
and capacitances. 

From the size of the Fabry-Perot diamond we have the 
following three equations, where we apply the source- 
drain voltage to the source electrode and keep the drain 
electrode at ground (see Fig.[TJc-f)). 



AE + 4Uc, 
AE + NUc. 



(1) 
(2) 

(3) 



Where e is the electron charge, AVg, ^^gi^ and AVsd 
determines the the size of the Fabry-Perot diamond as 
shown in Fig.jljc), and C^, Cg, Cd, and C are the ca- 
pacitance of the dot to gate, source, drain and the total 
capacitance. N is the equilibrium number of electrons 
added to the dot from a zero-bias resonance to the first 
positive-bias resonance, i.e., from position d to e in Fig.[l] 
N can be given in terms for of the tunnel barrier asym- 
metry {a = Ts/Td): 



r a + l 



(4) 



a can be found from the conductance at resonance: 
Go = 16a/(a + l)^e^//i. For the resonance indicated 



with an arrow in FigjTJb) we find Go ~ 2.6e^//i, a ~ 0.3, 
and TV - 3.1. 

From the width (full width at half maximum) at reso- 
nances in gate (Wg) and bias (Wgd) (see Figjljc)) we can 
set up the following two equations: 
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(6) 



where the second equation is a good approximation when 
the asymmetry of the capacitive or tunnel coupling is 
not too large (see appendix). N' is the number of elec- 
trons added to the dot between Vgate = ^^^/^ from 
resonance. We estimate N' by integrating a Lorentzian 
density of state for each energy level on the dot: 



N' 



Wg/2 



WJ2 
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^ {e + jAV,f+{\W,y 



de (7) 



where the sum should include an appropriate number of 
energy levels. If only one energy level is included {j = 0) 
N' = 2, but for increasing number energy levels included 
A'"' saturates at a higher number (since the tails of the 
other levels contribute). For the device analyzed in paper 
it saturates at N' ~ 2.5. 

By solving the equations above we can find expressions 
for the following parameters: 
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9 meV 
0.5 meV 
4.3 aF 

152 aF 
158 aF 
315 aF 
4.3 meV 
1 meV 
3.3 meV 



(8) 



We have in the right hand column estimated the parame- 
ters for the device analyzed in this paper [22 . Note that 
> F > Uc > Ao, where Ao ~ 0.11 meV is the super- 
conducting energy gap (see below). 

We now return to the measurements shown in Fig.jljb), 
where two parallel conductance ridges are observed at low 
bias due to the density of states in the superconducting 
electrodes. The separation between these two rides is 
4Ao/e, yielding Ao ~ 0.11 meV. In the following we fo- 
cus on measurements performed on and off zero-bias res- 
onance at the two indicated positions in Figjljb). Cur- 
rent versus junction voltage {IVj curves) off resonance 
for large scale voltages is shown in Fig.[2ja), where the 
black curve is with superconducting electrodes and the 
red curve is with a small magnetic field (150 mT) to sup- 
press the superconductivity. At high bias Vsd > 2Ao/e 
transport is governed by quasiparticle transport and one 
Andreev reflection processes yielding an excess current, 
while at sub-gap bias Vsd < 2Ao/e transport are gov- 
erned by Andreev reflections and supercurrent [4 . A 
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FIG. 3: (color online) Current versus junction voltage on reso- 
nance at position indicated in Fig.[T]^b) for three different tem- 
peratures. From left to right: 75 mK, 150 mK, and 300 mK. 
Four probe voltage bias measurement (circles), and fit (solid 
red hne) using Eq. ^ with Rj = 1 .IkVi, R = the tem- 

perature at which the curve is measured, and Ic — 4.8 nA, 
4.8 nA, and 4.6 nA from left to right. 



FIG. 2: (color online) Current versus junction voltage on (a- 
b) and off (c-d) resonance at positions indicated in Fig.[T]^b). 
(a) Black curve is with the electrodes in the superconducting 
state, and red curve with a small magnetic field (150 mT) ap- 
plied to suppress the superconductivity, (b) Close-up of the 
supercurrent branch from (a), measured with a voltage bias 
setup (circles) and a current bias setup (triangles), (c-d) De- 
pendence of the diffusive supercurrent branch on temperature 
(c) and magnetic field (d). 

close-up at very low bias voltages, shown in Fig.[2|b), 
reveals a pronounced supercurrent branch with finite re- 
sistance, a so-called diffusive supercurrent branch [13]. 
The black circles are measured with a voltage bias setup 
as shown in Fig.jlja), while the green triangles are mea- 
sured with a current bias setup (sweeping from negative 
to positive current). For voltage bias measurements we 
have observed no hysteresis or switching in the IV j curves 
at any gate voltages. But for current bias measurements 
switching and hysteresis are observed whenever the full 
/Vj-curve has local minima and maxima, as observed in 
Fig.|2|b). Such local minima and maxima will for current 
bias measurements lead to switching in voltage and result 
in a hysteretic /Vj-curve. To resolve the full /Vj-curve 
we have therefore used voltage bias measurements in this 
paper. 

In Fig.[2|c) and (d) we show the temperature and 
magnetic field dependence of the diffusive supercurrent 
branch, which we will analyze in the following. The zero 
bias slope of the diffusive supercurrent branch in Fig.|2] 
yields a resistance of the order kilo ohm. For a Joseph- 
son quantum dot junction with only two channels as for 
a nanotube the Josephson energy Ej = Mc/'^e can be 
comparable to the temperature of the cryostat. Ther- 
mal fluctuations will therefore lead to fluctuations in the 
phase difference across the junction, and consequently 
give a supercurrent branch with finite resistance. In or- 
der to dampen these phase fluctuations and thereby in- 
crease the size of the supercurrent branch, we have de- 
signed the environment of the superconductor-nanotube- 



superconductor junction as described in Ref. [13]. The 
quality factor for the junction is Q < 0.5, i.e., strongly 
damped. The full /Vj-curve for a damped Josephson 
junction including the external components (without Rj) 
was calculated by Ivanchenko and ZiPberman^lSj and 
used with great success by Steinbach et. al. [16]. Since 
this device has considerable current contribution from 
multiple Andreev reflections at sub-gap bias voltage we 
have to a rough approximation included a constant resis- 
tor Rj. The full /Vj-curve can then be calculated as^l3] 

rn/ ^ r t f h-r,i{Ej/kBT) \ Vj 
/(T^..) = /c/m( ^_^^^^^^^^^J +- (9) 

where In{x) is the modified Bessel function of complex or- 
der, and T] = {hVsd)/{2eRkBT). To plot I{Vsd) versus Vj 
instead of Vgd we can use that Vj = Vsd — RIiysd)' There 
are two fitting parameters in this theory, the tempera- 
ture dependent critical current Ic(T) and Rj. In Fig.js] 
we show three I versus Vj curves measured at the same 
gate voltage for increasing temperatures. From left to 
right: 75 mK, 150 mK, and 300 mK. The black circles are 
the measurement and the solid red curve is theoretical fit 
with Eq. (|9|. The three fits are made with Rj = 7.7kQ, 
and the temperature at which it is measured, the only 
free fitting parameter is Ic(T) yielding 4.8, 4.8, and 
4.6 nA respectively. Eq. ([9| fits the measured IV j curves 
very well for all temperatures with Ic as the only fitting 
parameter. Above ~ 300 mK smaller and smaller critical 
currents are needed to make a good fit. Critical cur- 
rents versus temperature found by these fits are plotted 
in Fig.[4]^a). At temperatures lower than ~ 300 mK the 
critical current is saturated at ~ 5 nA, while at higher 
temperatures it decreases more rapid than a BCS-gap 
dependence. In Fig. [4]^ a) we also plot the excess current 
versus temperature, measured at Vsd = 4Ao/e. We com- 
pare the measurement with theory for a superconducting 
quantum point contact [H |l9l |20l [21] . We use Eq. 1 and 
2 in Ref. [4 with A = A(T) having a BCS temperature 
dependence to fit the measured temperature dependence 
of the critical and excess current, solid red and blue curve 
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FIG. 4: (color online) (a) Temperature dependence of the 
measured critical current (squares) and excess current (dia- 
monds) on resonance, (b) Magnetic field dependence of the 
measured critical current on and off resonance (squares and 
circles), and excess current on resonance (diamonds). The 
normal state zero-bias conductance is 2.6e^/h on resonance 
and lAe^/h off resonance. The solid lines in both (a) and 
(b) are the predicted curves for a superconducting quantum 
point contact multiplied by a constant factor of 0.25 for the 
critical current and 0.7 for the excess current. Insert shows 
the magnetic field dependence of the sub gap structure of a 
similar device in the Coulomb blockade regime. 



m Fig.ga). The mamitude of the measured critical and 
excess current is 0.25 and 0.7 lower than the theory pre- 
dicts, while their qualitative dependence on temperature 
fits well with theory. 

In Fig.[4]^b) we plot the magnetic field dependence of 
the critical current on and off resonance (see arrows in 
Fig.jljb)), and excess current on resonance. The criti- 
cal currents are found by the same method as above by 
fitting Eq. ([9| to each measured IVj curve in Fig.[2|d). 
We compare the measurement to the same theory as 
above, but with A = A{B) = (1 - B/Bc)Ao, where 
Be ~ 90 mT is the critical field. We use a linear de- 
pendence because, as shown in the insert of Fig.[4]^b), 
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FIG. 5: (color online) (a) Differential conductance versus bias 
and gate voltage using Eq. (17) and (18). (b) and (c) Differ- 
ential conductance versus gate voltage at Vsd — mV (b), and 
versus bias voltage at resonance (c). Red squares are exper- 
imental data (measured with B — 150 mT), Solid black line 
is a Lorentzian fit to the measurement yielding Wg = 0.4 V 
and Wsd — 8-5 mV, and blue circles are numerical theory ex- 
tracted from (a). 



the sub-gap structure has approximately a linear depen- 
dence on magnetic field. The theory seems to fit qual- 
itatively well to the measurement. But the magnitude 
of the measured critical and excess current is, as above 
for the temperature dependence, 0.25 and 0.7 lower than 
theory. 

The dot is in the open regime with a charging energy 
oi Uc ^ 0.5 me V as discussed in the beginning of the 
paper, which is several times larger than the supercon- 
ducting energy gap (Aq ~ 0.11 meV). We speculate that 
the discrepancy of the factor 0.25 between the measured 
critical current and theory could be due to the charg- 
ing energy being larger than the gap thus suppressing 
the Cooper pair transport. The 0.7 discrepancy for the 
excess current has been seen before but we have no 
good explanation for that. 
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MEAN FIELD DESCRIPTION OF 
FABRY-PEROT RESONANCES IN A 
NANOTUBE QUANTUM DOT. 

The electronic states in the nanotube can be described 



by 



eff 



where 



"mrja 



f3=g,s,d 



c 



-UcN^-eVeffN (10) 



(11) 



to gate, source and drain. In the experiment we apply 
asymmetric bias, i.e., Vs = Vsd and Vd = 0. In the mean- 
field approximation (which is valid when F Uc), the 
Hamiltonian is 

Hcnt ^ ^ + UcN{N) - eV.ffN, (13) 

mrja 

where the total occupation (N) should be determined 
self-consistently 



(14) 



C 



(12) with 



and AE is the level spacing. The quantum numbers 
m, cr, T] describe the orbital, spin and pseudospin degrees 
of freedom, respectively. The subscripts g, s, and d refer 



E 



mrja/ / ^ p 

a=s,d 



27r 



nriu; + e 14)A^^^(cj). (15) 



J 



Assuming all levels to be simple Lorentzians with equal widths, the spectral functions are 



{uj-mAE- Uc{N) + e V^f^Y + {T /2f 



(16) 



Inserting this into the integral, summing over quantum numbers, and setting T = 0, then gives the self-consistency 
equation 



m a—s,d 



1 1 _J mAE + eVa + Uc{N) -eVeff 
tan I 



2 TT 



r/2 



(17) 



This is equation can be solved numerically. Once we know the total occupation for given gate, source and drain 
voltages, the current is given by 



e 4r,r, 



h r 



4e ATsTd 
h 7tT 



E 



tan" 



.1 / mAE + eVd + Uc{N) -eVeff 

[ r/2 



■ tan" 



.1 mAE + eVs + Uc{N) -eV.ff 

[ r/2 



(18) 



In Fig.jsja) we plot the differential conductance versus 
bias and gate voltage using Eq. (17) and (18) with the 
parameters found in Eq. (|8|. We compare the theory with 
experimental data measured with = 150 mT in Fig.jsjb) 
and (c). In (b) we make a gate trace at zero bias and in 
(b) we make a bias trace at the resonance indicated in 
Fig.Ba). 
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